Avalanche dynamics on a rough inclined plane 
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Avalanche behavior of gravitationally-forced granular layers on a rough inclined plane are investi- 
gated experimentally for different materials and for a variety of grain shapes ranging from spherical 
beads to highly anisotropic particles with dendritic shape. We measure the front velocity, area and 
the height of many avalanches and correlate the motion with the area and height. We also measure 
the avalanche profiles for several example cases. As the shape irregularity of the grains is increased, 
there is a dramatic qualitative change in avalanche properties. For rough non-spherical grains, 
avalanches are faster, bigger and overturning in the sense that individual particles have down-slope 
speeds u p that exceed the front speed Uf as compared with avalanches of spherical glass beads that 
are quantitatively slower, smaller and where particles always travel slower than the front speed. 
There is a linear increase of three quantities i) dimensionless avalanche height ii) ratio of particle to 
front speed and iii) the growth rate of avalanche speed with increasing avalanche size with increasing 
tan# r where r is the bulk angle of repose, or with increasing /3p, the slope of the depth averaged 
flow rule, where both r and f3p reflect the grain shape irregularity. These relations provide a tool 
for predicting important dynamical properties of avalanches as a function of grain shape irregularity. 
A relatively simple depth- averaged theoretical description captures some important elements of the 
avalanche motion, notably the existence of two regimes of this motion. 

PACS numbers: 45.70.Ht, 45.70.-n, 47.57.Gc 



I. INTRODUCTION 

Granular materials form phases with strong similari- 
ties with ordinary phases of matter, i.e., solids, liquids 
or gases In many natural processes, the coexistence 
of two of these phases is observed, requiring a complex 
multi-phase description. An example is avalanche for- 
mation, which occurs under various circumstances in na- 
ture (snow avalanches, sand avalanches on dunes, rock 
avalanches, land slides, etc.) as well as in industrial pro- 
cesses involving granular materials. Although laboratory 
realizations of granular avalanches do not have the full 
complexity of avalanches encountered in nature, labora- 
tory studies of this phenomena are important for under- 
standing the fundamental behavior of avalanches where 
there is a interesting combination of stick-slip friction, 
yield criteria for the solid phase, and the fluid-like mo- 
tion of the avalanche itself. Although the statistics of 
avalanche occurrences is a fascinating area of research 
Q , we focus here on the dynamics of individual avalanche 
events. 

There are several classes of experiments for avalanches, 
each with advantages and disadvantages. One is to slowly 
rotate a closed cylinder that is about 50% full of granular 
material 0, 0, IE 0]- As the cylinder rotates, the angle 
of the granular surface exceeds the critical angle # c , and 
material starts to flow along the surface. At rapid rota- 
tion rates the flow is continuous but for slower rotation 
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rates, avalanches occur because the rate of depletion of 
the granular material in the avalanche brings the surface 
back to the angle of repose 6 r faster than the rotation can 
maintain a surface angle greater than C . Significant ad- 
vantages of this approach are that there is no need to con- 
tinuously supply grains to the system, and the flow rate 
is easily controlled by the rotation rate. Further, using 
transparent boundaries enables direct observation of the 
velocity profile in the flowing layer because avalanches in 
this system typically extend to the side walls. From such 
optical measurements one finds that the vertical velocity 
profile of an avalanche follows an exponential decay |5[ 
in contrast to steady flows where it has an upper linear 
part. On the other hand, the role of friction at the side 
boundaries on the vertical velocity profile is unclear and 
may be quite important because the horizontal velocity 
profile is a plug flow with two exponential boundary lay- 
ers at the walls. 

A second realization of avalanches is to add granular 
material to a heap near its peak to induce granular mo- 
tion can be three dimensional or can 
be confined between rigid barriers, typically transparent 
glass or plastic plates, to form a unidirectional flow. For 
high input mass flux near the top, the grains flow continu- 
ously down the surface formed by other grains. For lower 
incoming mass flow, avalanches form intermittently. The 
transition between continuous and avalanching flow as a 
function of input mass flux is hysteretic 0, [Tcfl . Studies 
using diffusing- wave spectroscopy Q and molecular dy- 
namics simulations provided detailed properties of 
this intermittency. The experimental results show that 
the microscopic grain dynamics are similar in the con- 
tinuous and intermittent flow regimes [8]. In the con- 
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tinuous flow regime for a bulk heap, the grain velocity 
u(z) decreases linearly as a function of the depth z below 
the surface. Below some depth, however, a much slower 
"creep" motion is observed with an exponential velocity 
profile 0, El Bjj. For the case of avalanches, theoreti- 
cal results suggest that the avalanche amplitude should 
depend strongly on the velocity profile [14|, EH . 

A third laboratory system for avalanche studies is a 
layer of grains on a rough plane inclined at an angle 9 
with respect to horizontal 0, Ell, El, El, E3, El EI E3, 
E3, EB| • Because we use this system in the experiments 
reported below, we review some aspects of these inclined 
layer flows in more detail. Grains on a rough inclined 
plane form a thin stable static layer even for 9 > 9 C and 
the onset of the flow is expected only above a critical 
thickness h c , while the flow subsides at h s < h c [26L l27j|. 
The values of h c and h s are rapidly decreasing by in- 
creasing 9 as is schematically illustrated in Fig. [TJ Here 
h s is an average of the curves measured for 8 different 
materials (for the data see Section Hill) . and h c is an ap- 
proximate curve for illustration. 
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FIG. 1: (Color Online) The minimum layer thickness h s (av- 
eraged for 8 different materials) at which flow stops (solid 
line) and an illustrative curve for the critical thickness h c at 
which flow starts (dashed line). The operating conditions for 
avalanches (for the present work) fall in the region represented 
by the hatched zone. 

According to Bagnold [Hj], the shear stress in gran- 
ular flows is proportional to the square of the strain 
rate. This hypothesis has been checked in various con- 
figurations and has been proven to work relatively well 
for dense inclined plane flows [29]. The Bagnold stress- 
stain relationship leads to a velocity profile of the form 
u{z) = uq(1 — (h — z)/h) 3 / 2 where the surface velocity uq 
depends on h as uo ~ h 3 / 2 . The above convex velocity 
profile is recovered in MD simulations for relatively thick 
flows [30|, EH . As the internal velocity profile is difficult 
to determine experimentally, the most straightforward 
experimental test of the Bagnold hypothesis is to mea- 
sure the height dependence of the surface velocity uq or 
the depth- averaged velocity u. This relationship, called 
the flow rule (FR), was obtained for homogeneous flows 
with glass beads or sand [H, El, E3 an d is found to be 



consistent with the Bagnold flow profiles. Two slightly 
different forms of the FR are the Pouliquen flow rule 



(i) 



u/y/gh = l3 P h/h s (0) - 7 



and a modified form due to Jenkins (34| which we denote 
the Pouliquen- Jenkins flow rule: 



(2) 



u/y/gh = p PJ (h/h a (0))(tsm 2 9/ tan 2 #1) 



where 9\ is the vertical asymptote of the h s {9) curve. 

In the inclined layer system, one can either prepare the 
layer in a metastable state and mechanically induce a sin- 
gle avalanche [I?], El, EH E3 > or an intermittent series of 
avalanches can be produced by slowly adding grains near 
the top of the inclined layer as is done here. In the for- 
mer case the whole layer becomes metastable [13, EH, EH , 
and the shape and propagation of the avalanches depends 
critically on 59 (the level of metastability), giving rise to 
downward and also upward expanding avalanches when 
triggered with a small perturbation. For both kinds of 
avalanches on an inclined plane the angle at which flow 
starts or stops depends on the thickness of the layer con- 
trary to the heap experiments, on the properties of the 
bottom boundary such as surface roughness (which con- 
trols the zero- velocity boundary condition), and on the 
particular granular particle properties including inter- 
particle friction and particle shape. The distinct advan- 
tages of this system for the study of avalanches include 
the robust flow rule relationship in the continuous flow 
regime [26|, EH, EE[ that relates depth- averaged velocity 
and height, a simpler vertical velocity profile, and the 
ability to adjust the stability of the layer by changing 
the plane inclination angle 9. Note that both flow rules 
connect the rheology deep in the flowing state with prop- 
erties of the boundary between the flowing and the sta- 
tionary states. 

The majority of laboratory experiments use spherical 
beads - an idealized granular material. However, proper- 
ties of waves, either on a rough incline [33[ or in hopper 
flows (36[, depend on the shape anisotropy of the grains. 
Also, in recent numerical simulations the velocity pro- 
file in Couette flow depends strongly on the angularity 
of the particles [37[. Thus, it is natural to expect that 
avalanches will be similarly affected by the shape of the 
individual granular particles. Indeed, glass beads and 
sand show qualitatively different avalanching behavior 
[38( ] despite qualitatively similar flow rules in the steady 
flow phase. Avalanches formed by sand particles are 
larger with more dynamic grain motion than avalanches 
forme d by glass beads. By plugging the known flow prop- 
erties [33j into the depth-averaged model equations these 
basic differences can be explained (38[. 

In this paper, we present the results of an extensive 
study, using a set of different materials, and describe the 
details of our experimental methods used for the results 
presented in [38]. Our aim is to capture how the dy- 
namical properties of the avalanches depend on the grain 
shape irregularity, including the effects due to deviations 
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from the idealized spherical shape of the overall grain 
form and effects of the microscopic surface roughness of 
the particles. We find relations by which we can quanti- 
tatively predict major properties of avalanches as a func- 
tion of the angle of repose r or the slope of the flow 
rule (3p, both of which provide a measure of grain shape 
irregularity. 

We study the case where the layer is initially stable to 
small perturbations. As new grains are added to the layer 
at the top section (5%) of the plane, the granular layer 
becomes locally unstable in a cyclic manner leading to the 
intermittent formation of avalanches. These avalanches 
propagate down the plane on top of the stable static layer 
in a stationary manner, i.e., the shape and velocity of the 
avalanche does not change appreciably. The dynamical 
properties of these avalanches are studied far down the 
plane from where they are formed. We use eight differ- 
ent materials including rough irregular shaped particles 
of different sizes and spherical beads. The materials are 
characterized in Sec. [TT] where the experimental condi- 
tions are also described. In Sec. [1111 we discuss how the 
properties of the avalanches depend dramatically on the 
shape (spherical, or irregular) of the grains. We describe 
a simple theory of avalanche flow for glass beads in Sec. 
HVl Conclusions are drawn in Sec. [VI 



II. EXPERIMENTAL SECTION 

In this section we describe the experimental apparatus 
and the characteristics of the granular materials. The 
experimental techniques and apparatus have been de- 
scribed in detail elsewhere [38|, [39[ so only essential fea- 
tures are presented here. 



A. Experimental setup 

A sketch of the experimental setup is shown in Fig. 
[21 The granular material flows out of the hopper at a 
constant flow rate Q. The grains first hit a small metal 
plate that disperses the material so that the mass flux 
per unit width F = Q/W is relatively homogeneous in 
the y direction transverse (of total width W) to the in- 
clination direction x. The grains are deposited on the 
top 5% of the plane. The metal plate helps reduce pos- 
sible electric charging of the particles during the hopper 
flow. The granular layer reaches a critical state locally 
and an avalanche is formed and travels down the plane. 
The layer into which the avalanche moved has approxi- 
mate thickness h s . The stable static layer is prepared by 
releasing grains at the top of the plane and letting the 
system relax, or by letting the hopper run for a longer pe- 
riod of time. The glass plate has dimensions of 220 cm x 
40 cm with a rough surface prepared either by gluing one 
layer of the same particles onto it or by using sandpaper 
with different roughnesses. The two main control param- 
eters of the system are the plane inclination angle and 



the nature of the granular material used. The incoming 
flux and the roughness of the plane could be also varied. 
As we show in the following, within the limits presented 
below, variations in these latter parameters do not influ- 
ence the results presented in this paper. The dynamics 
of the avalanches are recorded with a fast video camera 
(up to 2000 frames per second) about 150 cm below the 
incoming flow (Camera 1). The vertical profiles of the 
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FIG. 2: (Color Online) Schematic illustration of the experi- 
mental setup. 

avalanches along their symmetry axes are recorded with 
Camera 2 using a vertical laser sheet to detect the height 
differences (Laser 1). With another laser sheet (Laser 2) 
slightly inclined with respect to the glass plate, 2D height 
profiles of avalanches are reconstructed. 

The incoming mass flow rate (per unit width along y) 
Fq is fairly constant for all measurements of one material 
and is in the range 0.02 < Fo < 0.2 g/cm-s depending on 
the material. The value of Fo is chosen by determining 
the flux F w needed to produce a wave state where the 
entire width of the layer is in motion. The properties of 
such waves were explored experimentally for glass beads 
[lH . The flux Fo is set to approximately F w /2 to obtain 
a reasonable number of distinct avalanches in the mea- 
surement area. For F < F w , adjusting Fo only affects 
the total number of avalanches and their size distribu- 
tion rather than the individual avalanche characteristics 
that are the focus of this paper. 



B. Characterization of the materials 

The shape of the granular particles has a large im- 
pact on the qualitative and quantitative behavior of 
avalanches in our system whereas other properties such as 
the mean particle diameter d have less effect. For exam- 
ple, despite some relatively low polydispersity for the ma- 
terials investigated (< ±30% except for fine sand which 
had ±50%), we do not observe any effects attributable 
to size-induced segregation. To show the nature of the 



shape anisotropy qualitatively, we first show in Fig. [3] 
images of the glass beads (d = 400 ± 100 /im), fine 
(d = 200 ± 100 iim) and coarse (d = 400 ± 100 /im) 
sand and salt particles (d = 400 ± 100 /im). The glass 



• > 









FIG. 3: Microscopic images of the materials: (a) sand with 
diameter d = 400 ±100 /im, (b) spherical glass beads with 
d = 500 ± 100 /im, (c) fine sand with d = 200 ± 100 /im and 
(d) salt with d = 400 ± 100 /im. 

beads are quite spherical although no attempt has been 
made to eliminate slightly aspherical particles, and no 
quantitative analysis of asphericity has been done. The 
sand has many irregular shapes ranging from approxi- 
mately spherical to very angular. Finally, the salt grains 
take on very angular shapes reflecting the cubic sym- 
metry of salt crystals. For the purposes of our studies, 
we assume that the grains do not change appreciably 
with time because the number of realizations of flows 
is limited, and the particle velocities are small. Thus, 
grains are not subjected to the repetitive collisions that 
might smooth their shape thereby changing the shape 
anisotropy or inter-particle friction. The copper mate- 

FIG. 4: Microscopic images of the copper particles of size 
d = 160 =b 50 jim and with packing fractions (a) r\ = 0.25, (b) 
j] = 0.33, (c) j] = 0.5 and (d) 77 = 0.63. 




rials are used to explicitly introduce a controlled shape 
distribution and explore the consequences of very differ- 
ent shapes on avalanche behavior. In Fig.JH we show im- 
ages of commercial (ECKA Granules GmbH & Co., KG.) 
copper particles (mean diameter d « 160 ± 50 /im) with 
very different average shapes. The particles range from 
highly irregular, dendritic-like shapes to almost spher- 
ical. The overall anisotropy is measured by the static 
packing fraction 77 for each shape, ranging from 0.25 for 
irregular particles to 0.63 for spherical copper particles. 
For comparison, the packing fraction for sand is about 
0.56 and for spherical glass beads is 0.63. 

The next thing to consider is how the different shapes 
and properties of the materials affect the basic granular 
flow properties, in particular 6 r and 6 C of the bulk ma- 
terial. We do so by observing avalanches on bulk heaps 
with a high speed video camera. As material is added at 
a very small rate to the top of the pile, avalanches form 
and propagate downwards intermittently. The basic con- 
figuration is shown in Fig. [5] where a bulk pile of sand is 
shown with a line indicating the determination of local C 
(just before an avalanche) or r (just after an avalanche, 
not shown). 




FIG. 5: Image of a sandpile - a heap of coarse sand particles. 
The scale is indicated in the figure and the solid line indi- 
cates the local slope used to determine C as shown or more 
generally r as well. 

The distribution of the critical angle C and that of the 
angle of repose r is plotted for 100 avalanches for each of 
the different materials investigated. The distributions for 
the glass, sand and salt are shown in Fig. [6] and the dis- 
tributions for the copper particles are presented in Fig. 
[3 The average values of C and r are indicated. As 
expected, 6 C and r are significantly higher for the piles 
formed by particles of irregular shape compared to piles 
of spherical particles. The widths of the distributions are 
also larger for irregular particles because shape irregular- 
ity gives rise to a larger variety of configurations and a 
wider range of angles at which avalanches start or stop. 

For thin layers on an inclined plane, grains start or 
stop flowing at angles determined by the layer thickness 
h as well as by their individual properties such as shape 
or surface roughness. In our experiments, we determine 
for each material the height h s (0) at which the material 
stops flowing, the thin layer equivalent of the bulk angle 
of repose. In a technique described in detail elsewhere 
(32| , we allow material to flow by slowly adding grains at 
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FIG. 6: (Color Online) Distributions of the critical angle C 
(filled columns) and the angle of repose r (open columns) 
measured for 100 avalanches on a three dimensional sandpile 
for (a) sand (red) and glass beads (black), (b) fine sand (red) 
and salt (black). 
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the top of the plane resulting in intermittent avalanches. 
Upon stopping the input, the layer comes to rest, and 
the volume of grains on the whole plane is determined. 
Knowing the surface area of the plate allows a determina- 
tion of the mean height h s with high precision, typically 
2%. Alternatively, a large quantity of grains is placed on 
the inclined plane, and a continuous flow persists over the 
entire plane for 10 - 20 seconds. After the flow subsides, 
the thickness of the resulting static layer is the same to 
within experimental error as that determined by the first 
method. 

Many of the details of the characteristics of flow on 
our inclined plane, including the detailed comparison 
with the Pouliquen flow rule [26[ and the effects of sur- 
face roughness, have been presented elsewhere [32]. We 
summarize our results for the particular materials used 
here by collapsing the data using a normalized ratio of 
tan#/ tan# r where r is determined from the bulk mea- 
surement and normalizing the height by the particle di- 
ameter d, see Fig. [U The continuous lines are best fits 
to the formula h s = ad/ (tan — tan#i). The resulting 
values for the fitting parameters a and 6\ are indicated 
in Table H For all the materials, the stopping height h s 
decreases with increasing as shown in Fig. dropping 
rapidly for angles close to the bulk angle of repose and 
more gradually for larger 0. The curves in Fig. [8] are 
for sandpaper with a roughness of R = 0.19 mm. The 
corresponding data curve for sand on a surface prepared 
by glued sand particles is identical within experimental 
error. The collapse of the data suggest that differences in 
plane roughness do not translate into significant changes 
in h s . We have also tested the influence of the incom- 
ing flux and kinetic energy of the incoming particles on 



FIG. 7: (Color Online) Distributions of the critical angle C 
(filled columns) and the angle of repose r (open columns) 
measured for 100 avalanches on a three dimensional sandpile 
for copper with packing fractions (a) 77 = 0.25, (b) 77 = 0.33, 
(c) 77 = 0.5 and (d) 77 = 0.63. Particle sizes d = 160 =b 50 jim. 



TABLE I: Fitting parameters a and #1, resulting as best fits to 
the data in Fig. [8] using the formula h s — ad/(tan0 — tan0i). 
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h s . For all the above configurations tested using different 
initial flow rates, the data points fall on the same curve 
within ±5% for tan 0/ tan r > 1.1. When approaching 
r the measurements become less accurate with the rapid 
increase of h s leading to a ±12% uncertainty of the data 
points for tan#/ tan# r < 1.1. 

There are several additional remarks about the system 
that are important. First, the majority of our measure- 
ments do not require detailed information about h c (0), 
the thickness at which grains start to flow for a particular 
0. Thus, we have not made measurements of h c . When 
data about this quantity are needed as in the modeling 
section presented below, h c (6) is estimated from mea- 
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FIG. 8: (Color Online) Static layer thickness normalized by 
d as a function of tan# normalized by tan# r . The continuous 
lines are best fits to the formula h s — ad/ (land — tan#i) with 
the resulting fitting parameters a and 0\ indicated in Table 
[U A vertical solid line indicates r . The horizontal bracket 
near the top of the figure indicates the range of over which 
avalanches are measured. 



surements in similar systems [26]. Second, avalanches 
are observed in the range of 1.04 < tan0/ tan# r < 1.34 
as indicated by the bar on the top of Fig. [HJ For higher 
0, even if a homogeneous static layer is prepared before- 
hand, the kinetic energy of the incoming grains is enough 
to slowly erode the pre-existing layer, and avalanches can 
only be observed for a short time. 



C. Experimental procedures and analysis 

For different materials and for different angles 0, we 
produce avalanches by adding grains continuously at the 
top of the plane. The main quantities of interest are the 
velocity of the front or of the individual grains, the lateral 
extent of the avalanche measured by the avalanche area 
A, and the height of the layer h as a function of position. 
These quantities are measured by analyzing images taken 
with the high-speed video cameras. Two images taken of 
transmitted light (Camera 1) for coarse sand at = 33.6° 
and = 38.1° are shown in Fig. [9l These avalanches are 
localized objects that can be isolated by differencing of 
subsequent images because regions were the grains move 
become visible on a uniform background. The area of the 
avalanche is determined by the areal extent over which 
finite displacements in image differencing are measured. 
Similarly the location of the front of the image differenc- 
ing region is used to determine the front velocity Uf. 

To track the motion of individual grains inside the 
avalanches, a sequence of images are extracted from 
movies of the granular flow dynamics [40]. Space-time 




FIG. 9: Images (a) at = 33.6° and (b) at = 38.1° using 
transmitted light taken with Camera 1 for sand with d = 400/x 
m. 



plots are created by taking the intensity along the sym- 
metry axis of the avalanche from movies taken by Camera 
1. An example is shown in Fig. [10] for sand. The vertical 
lines starting on the top of the diagrams are the traces 
of particles on the top of the static layer in front of the 
avalanche. The tilted trajectory dividing this region from 
the dynamic region is the trace of the avalanche front, 
providing a second method for determining Uf. The 
traces of moving individual grains inside the avalanche 
are tilted lines. The change in the tilt of these traces 
as we go downwards (towards the end of the avalanche) 
indicates the decreasing particle velocity u v behind the 
front. By analyzing the space-time plots the mean par- 
ticle velocity u p has been measured at several locations 
and the velocity profile u p (x) has been determined. For 
this example using sand particles, the particles in the 
avalanche core (close behind the front) move faster than 
the front velocity of the avalanche. 

In order to trace the height profile of avalanches, a ver- 
tical laser sheet (xz plane) is projected onto the plane. 
Movies are taken of avalanches that are cut by the laser 
sheet near their center (symmetry axis). The camera is 
mounted for these experiments on the side at an angle 
of about 13° with respect to the xy plane (Camera 2 in 
Fig. [2]). By measuring the speed of the avalanche and 
taking the intensity only along one vertical line of these 
movies, the profile of the avalanches can be traced assum- 
ing that the profile does not change as it passes through. 
In Figs. fTTla) and [TUb), profiles are shown for a sand 
and for a glass avalanche, respectively. The images are 
contracted by 25 times in the horizontal direction. This 
technique allows the instantaneous profile along one line 
to be determined. 

Another way to visualize the avalanches is to trace a 
laser line (Laser 2), see Fig. [2j across the direction of 
avalanche propagation, i.e., transverse to the plane inch- 
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FIG. 10: Space-time plot for a sand avalanche for — 36.8° 
and h s — 0.12 cm. Image size: 5.96 cm x 1.43 s. (Velocity 
profile shown in Fig. 120( a).) 



nation direction. Using the deflections of this laser line, 
the whole 2D surface of the avalanche is obtained. A sam- 
ple image (from Camera 2) is shown in Fig. [T2la) for a 
sand avalanche. From the image sequence, the 2D height 
profile is reconstructed and is shown in Fig. H2f b) for the 
same avalanche and similarly for an avalanche formed by 
glass beads in Fig. fT2tc). 



III. RESULTS 

The main objective of our experimental investigation 
is to understand the similarities and differences between 




FIG. 11: Avalanche profiles taken from the side with the help 
of a laser sheet for (a) sand (for = 33.6°) and (b) glass beads 
(for 6 = 22.6°). The horizontal size is contracted by a factor 
of 25. 



avalanches of granular materials with different properties, 
particularly grain shape irregularity. We first describe 
some of the qualitative characteristics of the avalanches 
and then provide quantitative measures of those proper- 
ties. 

The avalanches we investigated are stationary (do not 
change with time) above a certain size with an over- 
all shape that depends on the granular material. Glass 
bead avalanches have an oval shape, Fig. fl2lfc). whereas 
medium sized or large sand avalanches have two tails, 
Figs. and[T2lfb). which can break off and form very 
small avalanches. Small sand avalanches typically do not 
have these tails. Very small avalanches are not stationary, 
but decelerate, lose grains, and eventually come to rest. 
The length of the observation area is about 8-12 times 
the length of small avalanches. Avalanches that are de- 
celerating or that stop in the field of observation are not 
included in the data presented here. Interacting/merging 
avalanches are also eliminated. 

To demonstrate the stationarity of avalanche size and 
speed, the time evolution of the front velocity Uf and 
the characteristic size (the square root of the lateral area 
A) is shown for a set of avalanches in Figs. [T3l a) and 
fT3Tb) for sand taken at = 35.2°. The smallest, slow- 
est moving avalanches have a discernibly downward slope 
indicating a shrinking, decelerating avalanche. If the de- 
cay is exponential, the linear slope would yield a decay 
time of about 40 s for both the size and speed. For two 
avalanches of this initial size, one drops below a thresh- 
old, about 7 cm/s, and then quickly shrinks and decel- 
erates until it vanishes. Another with about the same 
initial conditions, decreases in size and speed for awhile 
but seems to recover and survive. 

To understand the difference in avalanche survival for 
small sand avalanches, one needs to take account of the 
response of the static layer to the passage of an avalanche. 



8 




FIG. 12: (Color Online) a) Image of the laser profile taken 
with Camera 1. The height of the avalanche is obtained as 
h — h s + 5x tan 0, where <j) is the angle between the plane and 
the laser sheet (Laser 2) and 5x is the displacement of the 
laser line. Height profiles of (b) sand avalanche for = 36.8°. 
Image size 6.9 cm x 34.6 cm, (vertical size rescaled by 25x) 
maximum height: h m = 0.35 cm, static layer thickness h s = 
0.12 cm; (c) glass bead avalanche for = 24.3°. Image size 
10.3 cm x 38.4 cm, (vertical size rescaled by 25x) maximum 
height: hm = 0.28 cm, static layer thickness h s = 0.18 cm. 



Medium or big sand avalanches typically take material 
with them along their center and deposit grains along 
their edges (see Figs. [9] and [T2]) . The change in the 
layer height is about ± 10% of h s . Thus, the different 
fates of austensibly similar sand avalanches suggests that 
small differences in residual layer thickness influence their 
evolution in that a small avalanche that encounters a 
slightly thinner region behind a recent large avalanche 
does not survive whereas a small avalanche propagating 
along the ridge formed by a large avalanche picks up a 
bit of mass and manages to survive (at least over the size 
of the interrogation window). No formation of residual 
ridges left by the wake of the avalanche is observed for 
small sand avalanches or for avalanches formed by glass 
beads. Although we have not studied this effect in detail 
for all the materials, one might surmise that the residual 
wake structure is a property of irregular grain avalanches 
whereas the spherical grain avalanches do not have this 
property. 

We now present a quantitative analysis of the speed 
and size of avalanches whose properties are stationary 
over the length of the channel, about 200 cm. The 
avalanche velocity Uf increases with avalanche size as 
shown in Figs. [TWa) and [TWb) for fine sand and glass 



15 



10 



a. 



-A A 
+ + 



AAA 
+ + + 



A 

+ + + 



V V V v V V 
t *> * © t ft 



o o 



3 



20 
15 
10 
5 



^ A A A A A 

+ + + + + + + + 



i 



V V 



a * I v ' v yv v v*w 



V 

* 

o o o 



o o 



_L 



1 2 3 4 5 6 

f(s) 

FIG. 13: (Color Online) Time evolution of the (a) lateral size 
A and (b) the front velocity Uf of four avalanches for sand 
with d = 400 /mi and = 35.2°. 



beads, respectively. Similar sets of data were obtained 
for all eight materials. For anisotropic grains, avalanches 
move faster with increasing see Fig. [T4t a). On the 
other hand, avalanches formed by spherical beads are in- 
dependent of in the sense that the curves obtained at 
various plane inclinations collapse, see Fig. [T4lf b). This 
observation is not in contradiction with data presented 
in [Hj], where avalanche speed for glass beads increases 
with decreasing 0. In fact, the tendency is the same here 
for the case of glass beads, as the average avalanche size 
increases with decreasing Fig. fl4lf b). so that bigger 
avalanches (with larger velocity) are observed for less 
steep inclines. A similar analysis for the anisotropic par- 
ticles is less conclusive, owing to the wide range of sizes 
and velocities measured. Nevertheless, after averaging 
Uf for all avalanches of different sizes, the Uf(0) curve 
for anisotropic grains is fairly independent of as well. 

To compare the Uf(A 1 ^ 2 ) curves for the different ma- 
terials, it is useful to employ dimensionless parameters. 
For gravity driven flows on an incline the relevant length 
scale is h s and the appropriate velocity scale is y/ gh s cos 
(4lj . The non-dimensionalization of the itf(A / 2 ) curves 
by these quantities collapses the data taken at various 
plane inclinations for each material. The collapsed curves 
are shown in Fig. [T5l a) for sand, salt, and glass beads 
and in Fig. [T5lfb) for the copper particles. The dimen- 
sionless avalanche velocity increases linearly with increas- 
ing dimensionless avalanche size for all materials. The 
data fall into two classes: irregular shaped particles form 
avalanches that reach sizes about two to four times larger 
than those formed by spherical beads. Further, the speed 
of the avalanches from the irregular grains is about 4 
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FIG. 14: (Color Online) Avalanche front velocity u/ as a 
function of the square root of the lateral avalanche area A for 
a) sand and b) glass beads. 
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FIG. 15: (Color Online) Dimensionless avalanche velocity as a 
function of avalanche size for (a) or sand (x), glass beads (o), 
salt (*), fine sand (A), and (b) copper particles with r\ — 0.25 
(A), rj — 0.33 (*), j] = 0.5 (x) and r\ = 0.63 (o). The dashed 
lines are best fits to the data, (c) Slope m u of the linear fits 
from (a) and (b) for all materials as a function of tan# r (*) 
or tan 0i (o). 



times larger than that of the avalanches formed by spher- 
ical beads. For a better characterization of these dif- 
ferences, we compare the slope m u of the best fits to 
the data (shown with dashed lines in Fig. [T5l a) and Fig. 
fTSTb)) for all materials. In Fig. H5f c) m u is plotted as a 
function of the tangent of the angle of repose r (or 6\ the 
vertical asymptote of the h s (6) curve), both of which pro- 
vide a measure of grain shape irregularity. We find that 
m u systematically increases with grain shape irregularity. 
In other words, the avalanche velocity increases with in- 
creasing avalanche size systematically for more irregular 
grains. Note that m u goes to zero at about tan# r = 0.36, 
implying that for materials with a small angle of repose 
(below about r = 19.8°) the dimensionless avalanche 
velocity is small (about Uf / ' \J gh s cosO w 0.35) and inde- 
pendent of avalanche size. 

Another measure of the avalanche size is its height, in 
particular the maximum height h m along the avalanche 
profile. We determine the height profile using the laser 
line method described in Sec. [Ill In Figs.fTTlfa) andfTTlfb), 
we see qualitatively that the height of glass avalanches is 
considerably smaller than that of sand avalanches. A sec- 
ond difference is that after the rapid increase in height at 
the front, a fast (exponential like) decrease is observed for 
sand avalanches whereas the height variation is straighter 



for avalanches formed by glass beads. 

Numerous avalanche profiles have been recorded for 
sand, glass beads and all copper samples at the same 
sequence of plane inclinations as before to provide a 
quantitative analysis of avalanche heights. Plotting the 
avalanche peak height h m /d as a function of the thick- 
ness of the underlying layer h s /d, we find a systematic 
linear increase as shown for glass beads and sand in Fig. 
[T6la), implying that at lower plane inclinations where the 
static layer is thicker the avalanche height is proportion- 
ally larger. The slope rrih of the curves gives a measure 
of the dimensionless avalanche height h m /h s for a given 
material which is basically independent of the avalanche 
size and plane inclination and is rrih = 1.45 and rrih = 2.5 
for glass beads and sand, respectively. Also shown is the 
slope m = 1.55 for glass bead waves on a velvet cloth [l8[. 
The close correspondence despite the considerable differ- 
ences in systems (cloth versus hard, rough surface) and 
phenomena (waves versus avalanches) is striking. For the 
copper particles (see Fig. fTBlfb)) rrih falls in the range of 
1.7 < rrih < 3.2 as r) changes between 0.25 < r) < 0.63. 
The general trend is that irregular particles form higher 
avalanches compared to avalanches of more spherical par- 
ticles. To quantify this trend we again plot h m /h s as a 
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FIG. 16: (Color Online) Avalanche peak height hm as a func- 
tion of the static layer thickness h s for (a) sand (x) and glass 
beads (o), (b) copper with 77 = 0.25 (A), 77 = 0.33 (*), 77 = 0.5 
(x) and 77 = 0.63 (o). (c) The ratio hm/h s as a function of 
tan# r (*) or tan#i (o). The continuous lines are linear fits. 




function of the tangent of the angle of repose r (or 0\ 
the asymptote of the h s (6) curve). As seen in Fig. fTBtc). 
h m /h s increases systematically with tan# r . This is an 
important relation as it provides a tool for predicting 
typical avalanche heights (for a given material) just by 
measuring the angle of repose of the material. 

Having demonstrated that both the dimensionless 
avalanche height h m /h s and the growth rate of avalanche 
speed with increasing avalanche size (defined as m u ) de- 
pend systematically on increasing grain shape irregular- 
ity, we can further explore whether such a tendency can 
be observed in other properties of avalanches. Using the 
space-time technique described in Sec. [Ill we can explore 
the ratio of the mean particle speed near the front u p to 
the front velocity Uf. 

To compare the behavior at the front, four examples 
are shown in Figs. H7Ta)-IT7lfd) for sand, copper with 
77 = 0.25, copper with 77 = 0.63, and glass beads, re- 



FIG. 17: Space-time plots taken at the symmetry axis of 
avalanches with camera 1. (a) sand at — 36.8° (9/6 c — 
1.01), (b) copper with 77 = 0.25 at = 39.5° (0/0 c = 1.01), 
(c) copper with 77 = 0.63 at 6 = 26.1° (0/0 c = 0.96) and (d) 
glass at = 23.5° (0/0 c = 0.95). 



spectively. In Figs. ITTT a) and ITTT b) one sees that for 
anisotropic particles the velocity of particles is larger 
than the front velocity; in many cases particles are flying 
out of the main body of the avalanche and are stopped by 
the static layer in front of the avalanche. On the contrary 
for the avalanches formed by spherical beads, Figs. flTTc) 
and fTTlf d) particles are slower than the avalanche front, 
and particles in the static layer begin moving just before 
the actual particles in the avalanche arrive at that posi- 
tion. This is clearly seen in Figs. fTTTc) andfTTTd) by the 
curved trajectories of particles near the front which are 
initially at rest. To quantify these qualitative differences 
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between particle and front motion, the average velocity of 
individual particles u p near the avalanche front and the 
front velocity Uf are measured for numerous avalanches. 
At a given plane inclination the particle velocity as 
a function of the front velocity shows a linear depen- 
dence, as indicated in the inset of Fig. [T8l The ratio 
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FIG. 18: (Color Online) Ratio of particle velocity u p and front 
velocity Uf as a function of the plane inclination for (a) sand 
(x) and glass beads (o), (b) copper particles with 77 = 0.25 
(A), 77 = 0.33 (*), 77 = 0.5 (x) and 77 = 0.63 (o). The inset 
of (a) shows Up as a function of Uf for sand (x) at = 36.8° 
and glass beads (o) at 6 = 25.2°. 

Up/uf is plotted in Figs. H5fa) and H5fb) as a function 
of the plane inclination for sand/glass-beads and copper 
particles, respectively. For glass beads and sand parti- 
cles, the separation is very clean with the spherical par- 
ticles moving slower than the front with 1 > u p /uf ~ 0.6 
whereas the irregular sand particles overtake the front 
with 1 < Up/uf ~ 1.3. There may be a slight downward 
trend with increasing for irregular particles but the 
data do not differentiate that trend from a near constant 
ratio. 

The ratios u p /uf are again less clearly separated for 
copper particles with ratios close to one for the 77 = 0.63 
and 77 = 0.5 particles. For whatever reason - per- 
haps inter-particle friction - spherical copper particles 
are marginal with respect to the separation of particle 
velocity and front velocity. The irregular copper parti- 
cles have Up/uf w 1.6, considerably greater than one, 
and are well separated from the more spherical particles. 



For the irregular particles there is a definite decrease in 
Up/uf with increasing 0. 

To quantify how the ratio u p /uf varies with grain 
shape irregularity, we again consider correlations between 
the ratio u p /uf and the angle of repose r (or #1). There 
is a linear increase of u p /uf as a function of tan# r (or 
tan 0i ) with a slope of about 3.5 as shown in Fig. [T9lfa) 
for the 6 materials presented in Fig. [181 We also compare 
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FIG. 19: (Color Online) Ratio of the particle velocity u p and 
front velocity Uf as a function of (a) tan# r (*) or tan#i (o) or 
(b) the slope (3 P j of PJFR Q or (c) the slope f3 P of PFR for 
sand, glass beads and copper with 77 = 0.33 and 77 = 0.5. In 
(b) the two data symbols correspond to the cases when f3pj 
was derived using r (*) or 6\ (o). 

how the ratio u p /uf varies with the equation relating u 
and h in the continuous flow regime, i.e., the granular 
flow rule. Earlier, we reported values of the slope (3pj 
in Eq. [2] for various materials [32] and found a system- 
atic increase of (3pj with increasing grain shape irregu- 
larity, i.e., with increasing 6 r (or 6\). Plotting u p /uf as 
a function of /3pj, we again find a systematic linear in- 
crease with a slope of about 2.8. The copper sample with 
77 = 0.5 seems to be anomalous, perhaps related to the 
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peculiar dynamics of this set of copper in that it is the 
only copper sample emitting strong sound during shear- 
ing. This emission is similar to but much stronger than 
the sand from the Kelso dune which is known to be an 
example of "booming sand dunes" [42| . Interestingly this 
peculiarity is not reflected in the dependence oiu p /uf on 
r , see Fig. HHa). 

Although the Pouliquen flow rule (Eq. [T]) does not per- 
fectly describe the flow properties of homogeneous flows, 
especially for the case of copper particles [32[ , for the clas- 
sification of avalanches we have determined /3p by taking 
only data for relatively shallow flows where h/h s < 10. 
Plotting Up/uf as a function of /3p (see Fig. fl9tc)) we 
again find a systematic linear increase. We will come 
back to this observation in Sec. IIV1 

Finally, we consider the detailed profiles of velocity and 
height along the propagating direction of the avalanche. 
The particle velocity u p (x) is obtained from space-time 
plots and the height profile from laser line measurements. 
The measurements are not made simultaneously but can 
be compared for avalanches with similar properties. We 
show profiles for glass beads and coarse sand in Figs. 
l20l a) and l20lf b) with scaled velocity and height on op- 
posite axes. The differences for sand and glass beads are 
striking. For glass beads, both velocity and height are 
approximately linear up to the maximum and then fall 
quickly over a steep but continuous front with a width of 
order 10h s . The slope of the surface behind the front is 
only slightly shallower than the unperturbed layer with 
SO/0 ~ —0.007 whereas the fractional angular increase 
near the front is about 0.1, i.e., a difference of about 3-4° 
relative to the plane inclination angle. 

For sand avalanches, a quite different picture emerges 
as shown in Fig.l20lfb). First, the velocity maximum and 
the height maximum occur at different values of down- 
stream distance with height peaking before velocity. Sec- 
ond, the height and velocity increase faster than linear 
from the back of the avalanche towards the front. Al- 
though we have adjusted the axes to align the maximal 
velocity and height of the avalanche, one can still see that 
the functional dependence for scaled velocity greater than 
one and h/h s > 2 are different. In other words, velocity 
and height are proportional for small values near the back 
of the avalanche but separate above values that, interest- 
ingly, differentiate in the mean between the spherical and 
irregular avalanche behavior. Despite the quite different 
width of the sand front, of order 25 h s , the fractional 
angular increase at the front is only slightly shallower 
than for glass beads, i.e., SO/0 w 0.07. The difference 
in profile for the sand arises from the nature of the par- 
ticles near the front. Whereas the glass beads form a 
compact avalanche with a sharp but distinct front, the 
sand particles overtake the front as defined by the mo- 
tion of the layer underneath creating a dynamics similar 
to a breaking wave. In the theory section below, we ex- 
plore a depth- averaged approach for explaining the differ- 
ent avalanche behavior of spherical and irregular granular 
particles. 
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FIG. 20: (Color Online) Avalanche height profiles (continuous 
lines) taken with the help of a laser sheet and velocity profiles 
(o) measured from space time plots for (a) glass beads (for 
= 24.3° and h s = 0.178cm) and (b) sand (for 6 = 36.8° and 
h s = 0.12cm). 



Although we cannot directly measure the interior pro- 
file of the avalanche, it is useful to provide a schematic 
illustration of the different type of avalanches to sum- 
marize what we have learned about them. In Fig. [2lJ 
we show the salient features of the two avalanche types. 
The case of shallower avalanches formed in materials with 
spherical particles is shown in Fig. l2Tta). Here the gran- 
ular layer fails, i.e., starts moving, just ahead of the 
avalanche (of order 5-10 h 3 ), leading to slower maximal 
particle velocities than the front velocity. The velocity 
and height of the avalanche are proportional and the 
variation of both quantities is linear in the down-plane 
direction. The form of the avalanche is reminiscent of 
a viscous Burgers shock as discussed in more detail be- 
low. In contrast, in Fig. l2JTb) a sand avalanche is visu- 
alized having a larger maximal particle velocity than the 
front velocity and a considerably higher avalanche height. 
The avalanche seems to behave like a "breaking wave" 
with low density particles spilling over the crest at speeds 
higher than the front speed. The material directly under 
the front seems to be solid with the avalanche more slowly 
entraining the underlying material. Direct measurements 
of this entrainment and the form of the avalanche pro- 
file would be very useful but are beyond the scope of the 
present work. 



IV. THEORY 

The theory of avalanche behavior is complicated by 
the existence of a moving fluid phase in contact with a 
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FIG. 21: (Color Online) Schematic view of avalanches formed 
in materials consisting of (a) spherical particles (e.g., glass 
beads) and (b) irregular shaped particles (such as sand). The 
horizontal scale (along the slope) is compressed by about a 
factor of 30 relative to the scale perpendicular to the plane. 
The tip region in b) represents low density material that spills 
over the front with speeds greater than the front speed. Note 
the difference in the mobilization of the underlying layer, es- 
pecially at the avalanche front. 



solid-like immobile granular phase. Thus, a complete the- 
oretical picture of avalanche behavior certainly requires 
a multi-phase approach. On the other hand, since the 
granular material is moving over most of its extent ex- 
cept for small regions at the front and back, a simpler 
approach neglects the solid state ahead and behind the 
avalanche. We can then use a depth- averaged approxima- 
tion, leading to the Saint- Venant shallow flow equations, 
adapted for granular media by Savage and Hutter [43j]. 
This approach was used in the analysis of wave forma- 
tion in dense granular flows [33[ and was applied by us 
to sand and glass-bead avalanches described previously 

0. 

For a flow of height h and mean velocity u, granu- 
lar flow down a plane, with the plane parallel to the 
x-direction, is described by the averaged equations for 
mass: 



dh 

7h 



djhu) 
dx 



and momentum balance: 
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= tanf 



H(u,h) 
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K — — gh cos 6. 

ox " 



The value of a is determined by the profile of the flow, 
a = 1 for plug flow (as in Ref. [43]), a = 4/3 for a 
linear flow profile, or 5/4 for a convex Bagnold profile 
[44j . The parameter K is determined by the ratio of the 
normal stresses in the flow: the stress parallel to the bed, 
cr xx , and that perpendicular to the bed, a zz . Numerical 
results show that K = <j xx /<j zz « 1 for steady-state flows 

The complicated part of this analysis concerns the fric- 
tion coefficient fi(u,h) [45[. For example, if the layer 
is too thin at a particular 0, the flow stops and a fric- 
tion coefficient appropriate for the flowing state is no 
longer valid. This transition from dynamic to static 
friction and the resulting yield condition for the solid 
phase in front of and behind the avalanche are not in- 
cluded in our approach and would be difficult to han- 
dle in a depth- averaged fashion. We will consider the 
consequences of this assumption later. For now we as- 
sume that the friction coefficient [i(u, h) is determined by 
the requirement that the steady flow obeys the rheology 
specified in Eq. ([5]), and will thus vary with the particle 
type. For simplicity, we consider only the Pouliquen flow 
rule; the Pouliquen- Jenkins flow rule adds algebraic com- 
plexity without improved understanding of this model of 
avalanche behavior. 

From the flow rule relation for steady-state flows (Eq. 
(PQ)) we obtain an expression for h s in terms of u and h: 



ha 



u + jvgh 



(5) 



The avalanches we will describe are observed for angles 
not very far from 6\ (i.e., tan#/tan#i < 1.4). In that 
range, the relationship 



h s 



ad 



tan# — tan 9\ 



(6) 



works very well in fitting the data for all the materials. 
The values for the two fitting parameters a and 0\ are 
shown in TableUfor sand, glass, and copper. As \i = tan# 
we obtain from (j6]) 



ad 

H — tan = tan 6\ + — . 

h* 



(7) 



(4) 



The slower avalanches we describe correspond to the 
case of low Froude numbers Fr = uj ^/ghcosO. For exam- 
ple, glass beads avalanches correspond to Fr w 0.4. Since 
the left hand side (LHS) of Eq. Q scales like Fr 2 w 0.16, 
we can set the LHS to zero. Plugging the expression for 
\i © into the right hand side (RHS), we get: 



ad T ^dh 

tan0 - tan 0i - K— = . (8) 

h s ox 

Substituting the expression for h s ((5|) into ([8]) gives an 
expression for u: 
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U= [tan 6 — tan 9\ — K^-) ■ (3 P h z ' 2 -^J~gh . (9) 
ad \ ox J 

Finally we use this form for u to substitute into (|3|): 



dh d 
dt dx 
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(10) 



which yields: 



<9i 9a; ad dx 
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(11) 



This equation has solutions similar to those of Burger's 
equation [46j . Thus, there is a solution consisting of a 
single hump propagating down the slope with velocity 
a(ft), with a smooth structure determined by the compe- 
tition between this nonlinear velocity term on the LHS 
and the dissipative term (RHS) of Eq. ([TT]) . Numerical 
solution of this equation yields a fairly linear ramp be- 
hind the front, very similar to Fig. l20la) for glass beads. 
Further, one can estimate the width of the front by bal- 
ancing the dispersive term and the dissipative term, tak- 
ing dh/dx ~ h/£ and d 2 h/dx 2 ~ h/£ 2 where £ is the 
front width. Inserting values for the parameters in Eq. 
([TT]) yields a scaled front width £/h s ~ 6, consistent 
with the height profile for glass beads in Fig. I20l a). We 
can also compare the predicted group velocity a (ft) with 
the measured avalanche front velocity for glass beads. 
Using the peak height ft m to evaluate a(ft), we obtain 
Uf/a(h) ~ 0.6, again consistent agreement between the- 
ory and experiment. 

The full equations of mass and momentum conserva- 
tion are hyperbolic with characteristic velocities 



a±y a(a — 1) 



K 



(12) 



If the velocity appearing in Eq. ([TT]) does not obey 
a < c+, then Eq. ([TT]) predicts a structure that moves 
faster than the maximum rate at which information can 
be propagated in the full system of equations, which is 
clearly impossible. For faster avalanches such as sand, 
the small Froude number limit is not valid since Fr > 1. 
In these circumstances, the Burger's type solution trans- 
forms itself into a truly discontinuous solution traveling 



at velocity c + [46], which is described by the full system 
Eqs. ([3E]) rather than by Eq. ([TT]) . 

We can directly test the conditions for which a < c+ 
is valid in terms of avalanche height. The normalized 
avalanche height h m /h s as taken from the slope of the 
fits in Figs. [T6l a) and [T6lf b) is shown in Fig. [22] for glass 
beads, sand and all four copper samples as a function 
of the flow rule slope /3p (using 7 = 0). We find a sys- 
tematic increase of h m /h s as a function of /3p, as we did 
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FIG. 22: (Color Online) Dimensionless avalanche height 
hm/hs (o) and the ratio of the particle and front velocities 
u p /uf (*) as a function of f3p for sand, glass beads and the 
four copper samples. The continuous line corresponds to the 
criterion a — c+, and the dotted line is a linear fit to the 
hm/h s (l3p) data. A horizontal dashed line denotes the value 

Of Up/Uf = 1. 

when plotting h m /h s as a function of tan# r (see Fig. 
fTBTc)). On the other hand, the boundary between the 
two regimes corresponding to a = c+ results in a critical 
thickness h / h c s riUcal of the flow that strongly decreases 
with increasing (3p (see the continuous line in Fig. [22]) . 
The range below h/h c s riUcal corresponds to smooth solu- 
tions of Eq. ([TT]) . with a < c+, whereas above h/h c s ritlcal 
theory predicts discontinuous solutions of the full system 
traveling at velocity c + . The line corresponding to the 
normalized avalanche height h m /h s as a function of (3p 
crosses this boundary, and the data points for sand and 
the two anisotropic (dendritic) copper grains fall above 
the curve, whereas data points for glass beads and the 
spherical copper samples fall below the boundary. This 
observation matches well with the result that the ra- 
tio of the particle velocity and the front velocity of the 
avalanche u p /uf also increases with increasing /3p and 
crosses the value 1 at about the same location corre- 
sponding to (3 P w 0.4 (see Fig. [22]) . 

Thus, we conclude that the glass bead avalanches, or 
more generally avalanches formed in materials with small 
/3p, reflect smooth solutions of Eq. ([TT]) . with a < c+, 
whereas avalanches in sand or copper samples with highly 
anisotropic grains, or more generally in materials with 
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higher (3p, represent discontinuous solutions of the full 
system traveling at velocity c+. The latter avalanches 
propagate into a quiescent bed because they are travel- 
ing at the characteristic velocity for the medium. In the 
experiments the signature of these "discontinuous" solu- 
tions is that grains move faster than the front speed, 
spilling or breaking ahead of the front in a low den- 
sity "foam-like" phase. The crossover between these two 
regimes happens at (3p « 0.4. The glass bead avalanches 
are analogous to "flood waves" in river flows, whereas the 
sand avalanches are analogous to "roll waves" in these 
flows [II, Ell. Note that ahead of the flowing avalanche, 
the moving material propagates into a material at rest, 
which is presumably in a state close to the critical Mohr- 
Coulomb state (48[. Unlike the flowing state, for which 
&xx ~ &ZZ-, in this critical state cr xx > <j zz . Thus the 
transition region, in which the flow accelerates from rest 
into a pseudo-steady state described by the continuum 
theory, can be viewed as a region of passive Rankine fail- 
ure, through which the compressive stress parallel to the 
bed, <j xx , is decreasing with time. The mechanics of this 
region are complex, and cannot be described by the Saint- 
Venant equations alone. 

We should point out that in the linear theory of the 
instability of steady flows, developed by Forterre and 
Pouliquen, the criterion a < c + corresponds to the stable 
regime of these flows with respect to wave disturbances 
[331 ] . Thus, our observation of discontinuous avalanches 
for sand and smooth avalanches for glass beads dovetails 
nicely with their observation that steady flows of sand are 
far more unstable to such disturbances than are steady 
flows of glass beads. 

The theory as presented has several nice features that 
semi-quantitatively describe the dynamics of avalanches 
in a range of materials with varying shapes and inter- 
particle friction. Nevertheless, the theory fails to account 
for several important aspects of the flow. First, the gen- 
eralized Burgers equation fails to capture the soliton-like 
nature of the individual avalanches, namely that they 
travel at constant speed and maintain uniform shape. 
One can include more physics by following Forterre and 
Pouliquen [45[ in providing a parameterization of the fric- 
tion coefficient that accounts for the static friction com- 
ponent and hysteresis between flowing and solid phases. 
Numerical simulations of either the full 2D equations or 
a ID version of the Saint- Venant equation with a more 
realistic friction, fail to resolve this difference [49j]. In 
particular, for glass beads one obtains a solitary localized 
solution but one with a flat rather than linear height pro- 
file behind the front. In that case, the speed is controlled 
by a front condition rather than by the mechanism con- 
sistent with our theory. We find empirically that the pre- 
dicted velocity a(h) is greater than the measured Uf, so 
perhaps a similar control is at play for the glass beads. 
For sand avalanches, full numerical simulations yield a 
form similar to our theory but again the height decays 
rather than assuming a constant profile. Thus, although 
the theory presented here captures some of the ingredi- 



ents of our experimental avalanches, the two phase na- 
ture of the avalanches is probably necessary for a more 
quantitative description of our results. 



V. CONCLUSIONS 

The main conclusion of our experimental study of 
avalanches is that the form of the avalanche depends on 
the granular particles. Differences are found not only 
at the quantitative level, but in the qualitative features 
(front propagation, avalanche shape, grain dynamics) 
which reflect the dramatic differences between smooth 
grains such as glass beads and irregular grains like sand. 
Measuring the basic avalanche properties such as veloc- 
ity and size, we find a systematic linear increase of the 
normalized avalanche peak height h m /h s with increas- 
ing tan# r or /3p, where r is the angle of repose of the 
material and (3p is the slope of the depth- averaged flow 
rule, both measures of grain shape irregularity. Simi- 
larly, the slope of the dimensionless velocity versus size 
curves increases linearly with increasing tan# r , suggest- 
ing that the avalanche velocity increases with increasing 
avalanche size systematically faster for more anisotropic 
grains. These two relations enable us to predict typical 
avalanche sizes and the growth rate of avalanche velocity 
as a function of increasing size for various materials, by 
simply measuring r or (3 p. 

When focusing on the dynamics of individual grains 
(and comparing it to the propagation of the avalanche), 
another characteristic difference is found for spherical 
and non-spherical particles. Particles with irregular 
shape move considerably faster than the avalanche front, 
indicating that inertia is an important ingredient in these 
faster moving avalanches. For avalanches formed by 
spherical particles, on the other hand, front propagation 
is transmitted through the contact points (force chains) 
near the front, thereby inducing motion through a col- 
lective yield (failure) condition at the front. In the lat- 
ter case, the velocity of individual grains is considerably 
smaller than the front velocity. When quantifying this 
important property of avalanches we again find a direct 
correspondence between u p /uf and either # r , the angle 
of repose of the material, or the slope of the depth aver- 
aged flow rule f3p (or /3pj). The value u p /uf increases 
linearly with increasing tan# r or /3p. These relationships 
again provide a tool for predicting the dynamical prop- 
erties of avalanches by simply measuring r or (3p for 
a given material. The direct correspondence should be 
valid for moderately polydisperse materials, but may not 
be valid for highly polydisperse materials where segrega- 
tion effects become important. 

In Section IIVI we have shown that a simple depth- 
averaged description can capture a strong change in the 
dynamical behavior of avalanches if we know the ba- 
sic (depth-averaged) rheology of the material. Also the 
striking difference between the behavior of avalanches 
formed by spherical and nonspherical particles warns us 
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that when modeling granular avalanches or granular dy- 
namics in general, predictions obtained by the simplest 
models neglecting shape irregularities might not be valid 
for systems with non- spherical particles. Our measure- 
ments, however, provide numerous additional quantita- 
tive features of avalanches for grains with different rhe- 
ologies that cannot be captured by the simple theory pre- 
sented here, leaving room for further theoretical or nu- 
merical studies that can provide a better understanding 
of avalanching behavior. 
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